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We propose a new observable, called zeal, to analyze events with jets in heavy ion collisions. The 
observable measures how a thermal medium affects the multiplicity and distribution of energetic 
particles in a jet. Using few known models for energy loss and jet quenching, we demonstrate its 
capability to distinguish the physics of these models. 


Since the seminal proposal [T] of using energy loss of 
fast particles and the related jet-quenching as a probe to 
study the nature of the hot medium formed in heavy ion 
collisions, a lot of experimental and theoretical studies 
have enriched our understanding of the medium produced 
in relativistic heavy ion collisions ElISl. For instance, 
Raa^ obtained from the ratios of single particle inclusive 
transverse momentum (p^)-spectra of nucleus-nucleus 
(AA) to that of suitably normalized proton-proton (pp) 
spectra shows a large suppression at RHIC HHS] and at 
the LHC [7]. Arguing it to stem from the leading par¬ 
ticle in the corresponding jets, this has been identified 
as the shining example of jet quenching by the medium, 
especially since no such suppression in seen in ratios con¬ 
structed for pA collisions. 

Several models of suppression for leading partons of a 
jet in a deconfined thermal medium have been proposed 
to account for these data. The essential picture is that 
the momentum of the leading parton transverse to its 
motion is broadened due to kicks from the medium. The 
broadening per unit distance {q = is related to the 
energy loss and is often used to characterize the interac¬ 
tion of leading partons with the thermal medium. 

The momentum broadening in the medium at a given 
temperature depends on the detailed description of the 
medium. Several groups [HHH] have calculated the sup¬ 
pression of leading partons in QGP . These involve dif¬ 
ferent models for the medium and the calculation of the 
energy loss. (For a comparison between the approaches 
see Refs. m-) Remarkably, a successful description of 
data seems feasible using models which treat the medium 
as a weakly coupled quark-gluon plasma (QGP), such as 
the GLV [8] formalism or the AMY formalism [12] as well 
as models based on the assumption that the medium is 
strongly coupled (sQGP), which treat q [9[ UHl E] as a 
non-perturbative parameter, or evaluate it using tech¬ 
niques suitable for strongly coupled theories [m. It 
is clearly desirable to find ways to distinguish between 
them. The extended pT-range at LHC is already a wel¬ 
come help in that direction. One may also look for new 
observables where these models, having already been con¬ 
strained by RaAi differ in their predictions for them. 

Energetic particles in quantum chromodynamics 
(QCD) come not isolated but as a part of a shower of 


several collimated particles known as jets. This is true 
already in the absence of the medium because of [HHni 
splitting of partons before fragmentation. In the QGP, 
interactions with the medium enhance these processes. 

LHC OUBHin] has ushered us in a new era of quantita¬ 
tive analyses of reconstructed jets, inviting thus a direct 
scrutiny of the original jet quenching ideas since one deals 
with truly fast particles. A jet is defined by clustering 
of observed particles in an event which depends on a pa¬ 
rameter that is correlated with the ang ular width o f jets 
or equivalently the cone radius R = y/ Ar]‘^ + The 

corresponding ratio R^2^2 depends, however, on the cone 
size (3- Intuitively this seems understandable mm 
since jets with larger R capture more momentum flow. 
However, it obscures a comparison with theory unlike 
Raa due to a i?-dependent variation in p^-suppression. 
Furthermore, one has to contend with a background sub¬ 
traction in heavy ion collisions, adding an extra source of 
difference in a comparison with pp collisions (or central 
vs peripheral collisions employed commonly). 

In this Letter, we propose a new observable for study¬ 
ing jets with the aim of minimizing, if not eliminating, the 
background subtraction and the cone size i?-dependence 
and yet capturing properties distinguishing the structure 
of jets in A A and pp collisions. We define transverse zea^ 
of a jet as follows, 

Z =-[log{Y^exp[-pT/{nT , (i) 

i 

where, the sum runs over the particles in the identified 
jet, pt is the magnitude of the total transverse momen¬ 
tum of the jet, fiT is the transverse direction of the jet, 
and ^ is the momentum of the ith particle in the jet. 
In the extreme case of only a single leading particle in 
the jet, Z = 1, whereas in the other extreme case of a 
jet of particles with typical energies of the order of the 
temperature T, Z 1. 


^ Here we use only the transverse momenta, hence the adjective 
“transverse”. Another definition, involving net magnitude of mo¬ 
menta, will give a different distribution. This may be useful 
to study the jet substructure involving energetic particles in pp 
collisions. Here on, we simply refer to Z as the zeal. 
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Since one expects that jets in A A collisions will typi¬ 
cally consist of a larger number of particles sharing the 
net pt due to broadening caused by both (a) medium 
induced bremsstrahlung and (b) thermal kicks in the 
medium, they will have a smaller zeal compared to the 
pp case. The difference may not be as stark as the ex¬ 
treme example above due to possible hadronization ef¬ 
fects. Nonetheless we expect that the zeal distribution 
observed in heavy-ion collisions to be shifted towards 
lower values compared to that in pp (or pA) collisions. 
Moreover, the shift will depend on the spectrum of parti¬ 
cles populating the jets in heavy-ion collisions compared 
to the spectrum in pp (or pA) collisions and will therefore 
encode quantitative information about how the medium 
affects the propagation of an energetic particle. As one 
needs a cone size R only to determine the pr of the jet, 
and as the background will be dominated by particles 
with much smaller transverse momenta, we expect the 
zeal distribution to be much less sensitive to both these 
factors. 

One may estimate the pt of the jet by using typical 
jet identification algorithms for the AA environment [23]. 
Alternatively, one can use the near side jet to identify the 
transverse momentum in the hard process and then an¬ 
alyze the zeal distribution of the away side jet. For a 
large enough cone radius i?, the value of the jet momen¬ 
tum pjet should not be particularly sensitive to R. The 
typical cone radius for pp collisions beyond which pjet 
saturates is about 0.6. R may be larger in A A collisions 
because of (a) momentum broadening transverse to the 
jet Pt due to the medium (b) greater probability of gluon 
emission. For a large enough i?, though, the sensitivity 
of the extracted value of pt should be reduced. 

Since the zeal weighs the contribution of multiple en¬ 
ergetic particles in the jet, it is more sensitive to the 
details of the distribution of the momentum between the 
constituent particles compared to the Raa which only 
follows the leading particle. Thus it certainly has a com¬ 
plimentary, perhaps greater, capability in differentiating 
models. For example, if a large fraction of the energy of 
the leading particle in an away side jet is lost to many 
thermal partons, the typical zeal values will be nearly 
zero. On the other hand if a few energetic partons carry 
away the energy then the typical zeal will still be finite. 
These scenarios may be relevant for the sQGP and wQGP 
respectively, subject to hadronization effects which them¬ 
selves will decrease the zeal. 

We test the viability of these ideas by investigating 
the difference between zeal distributions of pp and AA 
collisions for a few popular models of jet-quenching in the 
literature in different kinematic regimes. Our goal here is 


^ One appealing possibility is to use 7 or Z tagged jets, subject 
to the size of the initial state radiation m for the away side 
particle. 



FIG. 1. Distributions of partons in lOOGeV jets for vari¬ 
ous g as a function of The corresponding distributions 
for lOGeV jets are similar but narrower and peak at smaller 
See Ref. [25] for details. 

to elucidate the main features of the zeal distribution in 
heavy ion collisions and to check if it is a useful observable 
to characterize jets in the heavy-ion environment, and we 
find that it is. Computation of zeal of Eq. § needs the 
spectrum of particles which constitute a jet. While it 
can be taken from the experimental data or theoretical 
model simulations, we employ the published momentum 
distribution of particles in jets predicted by a few models. 

Let us first consider a model discussed in Ref. |25| . 
The paper describes implementation of a medium modi¬ 
fied QCD-splitting function in a Monte Carlo parton cas¬ 
cade for an energetic parton travelling through the fire¬ 
ball. The modification of the splitting depends on the 
jet-quenching parameter which is related to transverse 
broadening in the thermal medium. We chose it as a use¬ 
ful model since the spectrum of partons in a shower in 
both pp and AA collisions are published, allowing us to 
evaluate the zeal in the two cases and compare. We will 
also later consider some more popular models which are: 
GLV (based on [H [26l [27]), ASW (based on [23|29]), 
AMY (based on [9l l30l [3T] b 

One can obtain the detailed spectrum of particles in a 
jet in this model by simulating events using the Monte 
Carlo QHERWIG and identifying jets, as in [25]. Eor 
simplicity, here we generate an ensemble of jets as fol¬ 
lows. Consider a parton shower with total energy E. In 
QHERWIG, this is created by injecting a gluon with en¬ 
ergy E into the medium. We partition this energy into 
partons carrying a fraction x = exp(^) chosen randomly 
according to a probability distribution function which is 
identical to the predicted distribution function of partons 
given in Ref. [25] (see Eig. 4 in Ref. [25]). This ensures 
that to a good approximation the average (over a large 
number of jets constructed in this manner) distribution 
of partons as a function of the fractional energy carried is 
the same as the predicted distribution using QHERWIG. 
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FIG. 2. Zeal distributions for lOGeV and lOOGeV jets for various q. 


These distributions are given in Ref. [25] for E = 
lOGeV and E = lOOGeV and for qL = 0, gL = 2GeV^, 
qL = 5GeV^, qL = 20GeV^, qL = 50GeV^, where L is 
the path length of the parton in the medium. Using the 
software “gSdata” we extracted the data from the plots 
and employ as the probability density function for the 
procedure described above. Our resultant distributions 
of partons in a jet of energy lOOGeV, averaged over an 
ensemble of 200, 000 parton showers is given in Fig.[^ A 
comparison with Fig. 4 in Ref. [25| shows that the ob¬ 
tained distribution is similar to the predicted distribution 
in QHERWIG. The same is also true for the distributions 
for lOGeV jet energy which we have not shown here. 

If we consider events of equal centrality (or equivalently 
system size) one can study the range of medium-jet cou¬ 
pling from weak to the strong coupling. Alternately, if 
we know the value of g, the range of simulations covers 
central to peripheral events (assuming that the effect of 
the expected reduction in g due to a slight lowering of the 
temperature for more peripheral collisions is not substan¬ 
tial and the dominant effect is the change of geometry). 

Fig-i displays zeal-histograms of events for a single 
gluon jet in the medium with increasing g for both 10 
and lOOGeV. The case with g = 0 corresponds to the pp 
collision where the distribution for ^ is sharply peaked 
at 0.765. This is associated with jets where two particles 
share the bulk of the energy of the jet. Likewise, there 
is also a three-particle multiplicity peak clearly visible in 
the data at about Z ^ 0.54 in particular for lOGeV. We 
expect these few particle structures to be smoothened 
out in detailed simulations when hadronization is taken 
into account. 

We note that in Fig. smaller zeal values get popu¬ 
lated for both lOGeV and lOOGeV jets, as the medium 
broadening effect increases. In addition, the height of the 
two-multiplicity peak at Z = 0.765 decreases. The differ¬ 
ence between gL = 0, gL = 20 and qL = 50 scenarios is 
easily noticeable. The broad peak shifts towards smaller 
zeal values becoming somewhat narrower at larger qL 




FIG. 3. Ratios of distributions for lOGeV (upper panel) and 
lOOGeV (lower panel) jets with the distribution for gL = 0 
distribution for various g. 

values and the two-multiplicity peak vanishes. We ex¬ 
pect the main feature that the zeal distribution moves 
towards smaller values with larger qL to be robust and 
worth exploring further. 

In order to contrast the qL = 50 and gL = 0 cases bet¬ 
ter we construct their ratios, as shown in Fig. These 
may be recognised to be similar to the Rcp^ where one 
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studies centrality dependence of the medium effects, or 
Raa which are main tools to study jet-quenching exper¬ 
imentally since systematic effects cancel out in them. As 
for Rcp, no change in zeal distribution due to medium ef¬ 
fects will correspond to a constant line at unity for the ra¬ 
tio. We observe a net decrease in events with zeal values 
above Z - 0.4 {Z - 0.5) for E = lOOGeV {E = lOGeV) 
and net increase in zeals below, with increasing values of 
qL. For a fixed g, this amounts to increasing the path 
length or more central collisions. We also note that the 
ratio decreases as the zeal increases. In particular, very 
interestingly, the ratio of the few particle structures fall 
with increasing qL in a manner similar to the smooth 
region of the plot. 

Different models of medium effects will have different 
predictions for the zeal yields. Comparison with exper¬ 
imental data for the zeal distribution can possibly rule 
some out or verify them. It is well known (for eg. [13]) 
that the value of q may not be a good way to charac¬ 
terize different models because the same value of q gives 
rise to different medium induced bremsstrahlung rates in 
different formalisms, and consequently different values of 
Raa- It may be better to use the zeal distributions for a 
fixed centrality to differentiate various theoretical models 
of jet quenching. The medium induced bremsstrahlung 
spectra are directly related to those of relatively ener¬ 
getic particles moving along the jet direction, at least 
before hadronization. Therefore, it may turn out better 
to compare the zeal distributions predicted for various 
models for the respective medium effects tuned to obtain 
the correct Raa and see whether they can further assist 
in distinguishing the models by a comparison with the 
experimental results. 

Fig. 19 of Ref. [13] provided us the necessary input for 
our procedure above. It compares the models [i|9l|26]- 
[3T] by tuning the model parameters such that the value 
of R 7 — the suppression of a pr spectrum falling with the 
power 1/p^ — is 0.25 when the leading parton traverses a 
QGP brick of length L = 5fm. We generate ensembles of 
events for the initial injected energy of 20GeV partitioned 
into partons carrying a fraction x chosen according to the 
distribution given in Ref. [13]. Fig. [^exhibits our results 
for the zeal distribution for different models obtained by 
following the same procedure as described above. We 
find it very satisfying that even when the models are 
tuned to give the same suppression, the zeal distribution 
is able to distinguish between them, in particular showing 
a significant difference between the AMY and the GLV 
models. 

In summary, we propose a new jet observable, the zeal, 
defined by Eq. 0 to unravel the medium effects. It de¬ 
pends on the transverse momentum distribution of par¬ 
ticles in the jets, and is therefore more discriminating 
than Raa of the leading partons. Its advantage is that 
it weighs the energetic partons more heavily and hence 
is particularly sensitive to the processes that lead to the 




FIG. 4. The upper panel shows zeal distributions for 20GeV 
for various models described in Ref. m- The lower panel 
shows ratios of the other models with AMY. 

energy loss of the leading partons. For frequent medium 
induced bremsstrahlung with several gluons carrying a 
tiny fraction of the energy of the leading partons, the 
peak of the zeal distribution should move towards lower 
zeal values unlike the case where induced bremsstrahlung 
is rare and the emitted gluons carry significant fractions 
of the energy of the leading partons. 

Its second advantage is likely to be its smaller sensi¬ 
tivity to the background and the cone radius R. Since 
one can use large values of jet cone radii R to extract 
the pt of the jet to define zeal, it may be further useful 
since jets are expected to be wider in A A compared to 
pp collisions due to broadening. This will reduce the sys- 
tematics associated with the extraction of the pr of the 
jets, and hence the calculation of zeal. It is also infrared 
safe (an extra soft parton i will not change the zeal sig¬ 
nificantly because its contribution goes as exp(—p^/p^)) 
and is not affected significantly by accidental inclusion of 
some particles which are part of the thermal medium. 

Our results suggest that the zeal distribution is sen¬ 
sitive to medium effects, shifting it to smaller zeal val¬ 
ues with increasing centrality. We also see that it can 
be used to distinguish models of jet quenching. It will 
be interesting to study how these results are affected by 
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hadronization. We are currently pursuing this study and 
its application to the experimental data. 
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